Logical XOR gate response in a quantum interferometer: A spin dependent transport 
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We examine spin dependent transport in a quantum interferometer composed of magnetic atomic 
sites based on transfer matrix formalism. The interferometer, threaded by a magnetic flux <j>, is 
symmetrically attached to two semi-infinite one-dimensional (ID) non-magnetic electrodes, namely, 
source and drain. A simple tight-binding model is used to describe the bridge system, and, here 
we address numerically the conductance-energy and current-voltage characteristics as functions of 
the interferometer-to-electrode coupling strength, magnetic flux and the orientation of local the 
magnetic moments associated with each atomic site. Quite interestingly it is observed that, for 
<f> = <fro/2 (<f>o = ch/e, the elementary flux-quantum) a logical XOR gate like response is observed, 
depending on the orientation of the local magnetic moments associated with the magnetic atoms 
in the upper and lower arms of the interferometer, and it can be changed by an externally applied 
gate magnetic field. This aspect may be utilized in designing a spin based electronic logic gate. 
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I. INTRODUCTION 

With the rapid advancement in nanoscience and nan- 
otechnology, specially in nanofabrication techniques [l[, 
study of spin dependent transport in mesoscopic Q sys- 
tems has emerged as one of the most challenging top- 
ics in the last few decades. Analysis of spin transport 
and spin dynamics is essential to understand and develop 
the field - 'spintronics' @, |3]- With the discovery of Gi- 
ant Magneto-resistance (GMR) based magnetic field sen- 
sors Q in 1994, remarkable development has taken place 
in the field of magnetic data storage applications and 
quantum computation techniques. A drastic enhance- 
ment in computation time has been made possible using 
the idea of quantum coherence and spin entanglement. 
Manifestation of coherence is one of the most important 
aspect of mesoscopic systems. It is evident from theoreti- 
cal f6T-[l2j and experimental [l3l - fl5| studies of spin trans- 
port through quantum confined nanostructures that the 
conductance of such systems depends on the spin state 
of electrons passing through the system and it can be 
controlled by an externally applied magnetic field. But 
measurement of current through these ID nanostructures 
does not reveal the feature of quantum coherence, as it 
is detectable through interference experiments, most no- 
tably Aharonov-Bohm (AB) interferometry fll|. In order 
to study the effect of coherence, spin dependent transport 
has been studied in various types of ring type conduc- 
tors or two path devices 17( such as an AB ring or AB 



type interferometer with embedded quantum dots, with a 
magnetic flux (f> penetrating the area enclosed yielding a 
flux dependent spin transmission probability. The study 
of spin dependent transport through intcrfcromctric ge- 
ometries are important for further development in quan- 
tum information processing as well as for designing spin 
based nano-devices. The key idea of designing spin de- 
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pendent nano-electronic devices is based on the concept 
of quantum interference effect [l8l - [23j . and it is gener- 
ally preserved throughout the sample having dimension 
smaller or comparable to the phase coherence length. In 
realistic situation, experimentally sizable rings are typi- 
cally of the order of 0.4-0.6 fim. Therefore, ring type con- 
ductors or two path devices are ideal candidates where 
the effect of quantum interference can be exploited [24| . 

Recently, spin transport through AB type interferom- 
eters with embedded quantum dots has drawn much at- 
tention because of its demonstration of several physi- 
cal phenomena e.g., quantum phase transitions, resonant 
tunneling and many body correlation effects. It opens a 
new area of study of spin transport, which includes spin 
dependent conductance modulation, spin filtering, spin 
switching, spin detecting mechanisms, etc. Conductance 
of such mesoscopic systems is associated with the trans- 
mission probability (T) of electrons, which can be calcu- 
lated numerically by several methods like, mode match- 
ing techniques [25| , Green's function approach (2fj| - [28| or 
transfer matrix method [2!| [3(| • 

Aim of the present paper is to study the spin depen- 
dent transport through an AB type intcrfcromctric ge- 
ometry made up of magnetic atomic sites. The inter- 
ferometer, threaded by a magnetic flux cf>, is attached 
symmetrically to two ID semi-infinite non-magnetic elec- 
trodes. A simple tight-binding Hamiltonian is used to 
describe the system where all the calculations are done 
using transfer matrix formalism. Spin dependent conduc- 
tance is calculated using the Landauer formula and also 
current- voltage characteristics are computed through the 
Landauer-Buttiker formalism [3ll . [3^ |. We explore several 
important features of spin transport with this simple, yet 
interesting geometry. Quite nicely we see that, at the half 
flux-quantum value of (f> ((f) = (f>o/2), the system exhibits 
XOR gate like response depending on the orientations of 
the local magnetic moments in the upper and lower arms 
of the interferometer, that can be changed by an external 
magnetic field. To the best of our knowledge, the spin 
based XOR gate response in such a simple geometry has 
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not been addressed earlier in the literature. 

We organize the paper in this way. Following the intro- 
duction (Section I) where we address some general fea- 
tures and recent theoretical as well as experimental stud- 
ies on spin transport, in Section II, we describe the model 
and theoretical formulations for the calculation. In the 
theoretical formulation we describe in details the trans- 
fer matrix method following the renormalization proce- 
dure [H, [34[. Section III explores the numerical results, 
where we show the variation of conductance as functions 
of magnetic flux (</5), orientation of the local moments 
embedded in the intcrfcromctric arms and energy of the 
injecting electrons and then we illustrate the current- 
voltage (I-V) characteristics that clearly signify the log- 
ical XOR gate response. Finally, we summarize our re- 
sults in Section IV. 



II. MODEL AND SYNOPSIS OF THE 
THEORETICAL BACKGROUND 

Let us begin by referring to Fig. [TJ where a quantum 
interferometer, penetrated by an AB flux (f> (measured 
in unit of the elementary flux-quantum 4>q = ch/e), is 
attached symmetrically to two non-magnetic electrodes, 
viz, source and drain. The magnetic conductor i.e., the 
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FIG. 1: (Color online). Schematic view of a quantum inter- 
ferometer, threaded by an AB flux (j>, attached to two semi- 
infinite ID non-magnetic electrodes, viz, source and drain. 
The filled purple and blue circles correspond to the magnetic 
and non-magnetic atoms, respectively. 

interferometer is composed of six magnetic atoms, four 
of which are placed at four different corners of the AB 
type interferometer and the rest two are connected to the 
electrodes directly. The side-attached electrodes consist 
of infinite number of non-magnetic sites labeled as 0, — 1, 
—2, . . ., -co for the left electrode (source) and 7, 8, 9, . . ., 
oo for the right electrode (drain). Each magnetic atomic 
site has a local magnetic moment associated with it. The 
direction of magnetization in each magnetic site is cho- 
sen to be arbitrary and specified by angles 9 n and <p n 
in spherical polar co-ordinate system for the nth atomic 
site. Here, 9 n represents the angle between the direction 
of magnetization and the chosen Z axis, and ip n repre- 
sents the azimuthal angle made by the projection of the 
local moment on X-Y plane with the X axis. 

The Hamiltonian for the full system i.e., the electrode- 



interferometer-electrode can be described as, 

H = H D +H L +H R + H LD + H DR (1) 

where, Hp corresponds to the Hamiltonian of the AB 
type interfcrometric device made up of magnetic atomic 
sites. -ffL(i?) represents the Hamiltonian for the left 
electrode i.e., source (right electrode i.e., drain), and 
Hld(dr) is t ne Hamiltonian representing the device- 
electrode coupling. 

The spin polarized tight-binding Hamiltonian for the 
interferometer can be written within the non-interacting 
electron picture in the form, 

6 5 

H D = ^2 c n ( £ o - h n .<?) c n +^ (cjte 10 c i+ i 

n=l i=1 

+ cJ +1 te- i0 Ci) + (c{tc 2 + 4tCi) + 
(4tc 6 + cj,tc 4 ) (2) 

where, 
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In Eq. @, 1st term corresponds to the effective 
on-site energies of the interferometer. eo' s are the site 
energies, while the h n .cr term represents the interaction 
of the spin (a) of the injected electron with the local 
magnetic moment placed at the site n with strength h n . 
9 n and tp n represent the orientation of the local magnetic 
moment situated at the site n as mentioned earlier. 
This term is responsible for spin flip scattering at the 
sites. Flipping of spin violates spin conservation in the 
transport process through the magnetic conductor which 
may provide much impact in spintronic applications and 
we will discuss about it in the forthcoming sub-sections. 
Second term describes the nearest-neighbor hopping 
integral between the sites, at the corners of the interfer- 
ometer, modified due to the presence of AB flux 4> which 
is incorporated by the term O = 27r0/4<73o. The 3rd 
and 4th terms represent the nearest-neighbor hopping 
between the atomic sites 1, 2 and 4, 6, respectively. 
Similarly, the Hamiltonian -Hl(.r) can be expressed as, 

H L(R) = c i £ L(R-) c i + (^LWCi+l + h - c -) ( 3 ) 

i i 

where e_L(i?)'s are the site energies of the electrodes 
and ^l(_r) is the hopping strength between the nearest- 
neighbor sites of the left (right) electrode. In this Hamil- 
tonian, £l(R) an d tL(R) are in the form, 

( ZL(R) \ 
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In the same fashion, the conductor-electrode cou- 
pling Hamiltonian is described by, 

HlD(DR) = ( c o( 6 )tLD(DR) c l(7) + c i( 7 )tLD(DR) c 0(6) 

(4) 

where, t LD ^ DR ^ being the device-to-electrode coupling 
strength. 

Now, we start with the Schrodinger equation, 



where, 



l*) = S)[^t|it> + V'dH)] 



(5) 



(6) 



Here, |$) is expressed as a linear combination of spin up 
and spin down Wannier states. 

In order to calculate the spin dependent transmission 
probabilities through the interferometer, first we map 
the two-dimensional (2D) geometry into ID structure by 
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FIG. 2: (Color online). Schematic diagram of the renormal- 
ized ID geometry. In this renormalized geometry the total 
number of atomic sites N = 4, and, here we label the atomic 
sites in such a way (viz, 1, 2, 4 and 6) to understand the 
renormalized version of Fig. [T] much clearly. It shows that the 
atomic sites 3 and 5 of Fig. [1] gets renormalized. 

renormalization procedure. The schematic view of the 
ID geometry is shown in Fig. [2j 

We start renormalizing the intcrferomctric geometry 
by writing down the difference equations at the six sites 
of the interferometer. They are given as follows, 
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Substituting ips and ip5 using Eqs. (|9|) and (|11[) we obtain 
the renormalized difference equations for sites 2 and 4 as, 



(E - e 2 )V> 2 = tVi + t' 24 V>4 
(E - e 4 )V> 4 = ti/j 6 + t 4 2^2 



(13) 
(14) 



Due to renormalization, the site energies of the 2nd and 
4th sites get modified, and they are, 

£ 2 = £2 + t23 • (E — e 3 ) 1 • t 32 

+ t 2 5-(E-e 5 r 1 -t52 (15) 

£4 = £4 + t43 ' (E - £3) 1 • t 34 

+ t 4 5-(E-e 5 ) _1 -t 54 (16) 

The hopping term between these sites (2 and 4) is also 
modified as, 



t 7 2 4 — t 2 3 • (E — £3) • t3 4 

+ t 25 • (E - eg)" 1 ' t 54 



(17) 



and t 42 is the hermitian conjugate of t' 24 . 

With this renormalized ID geometry, we use transfer 
matrix method to calculate spin dependent transmission 
probabilities (T) and current-voltage (I-V) characteris- 
tics through the bridge system. 

For an arbitrary site n, the transfer matrix (P) can be 
defined in terms of the wave amplitudes of its neighboring 
(n + 1) and (n — 1) sites as, 




= P 




(18) 



In our case, the transfer matrix equation for the renor- 
malized geometry relating the wave amplitudes at sites 
0, -1 and N + 1, N + 2 becomes, 



V lpN+1,1 




(19) 



where, N corresponds to the total number of sites in the 
ID magnetic device after renormalization process. For 
our renormalized ID system, N = 4. M being the trans- 
fer matrix for the full system and it can be expressed 
as, 



M = M R -P 6 -P 4 -P 2 -P 1 -M L 



(20) 



where, Pi, P2, Pa and Pq represent the transfer matrices 
for the sites labeled as 1, 2, 4 and 6, respectively. Ml and 
Mr correspond to the transfer matrices for the boundary 
sites at the left and right electrodes, respectively. 



4 



To evaluate the transmission probabilities of up and 
down spin electrons, we calculate the explicit form of M, 
determining all the transfer matrices (Mr, Pq, P 4 , P2, 
Pi and Ml). The matrices are given below. 
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where, 



Oil 

a 
Pi 

$2 



4 4 [2,2](£-e 2 [l,l])+4Jl,2]e 2 [2,l] 

4[l,2]t^[2,2]+t / 24 [l,2](E-e / 2 [2,2]) 

t 24 [l,l]t 24 [2,2]-t 24 [l,2]t 24 [2,l] 

e 2 [2,l]t 24 [l,l]+4 4 [2,l](i?- e2 [l,l]) 

i 24 [l,l](£-e 2 [2,2])+^ 4 [2,l]e 2 [l,2] 

i 24 [2,2]4 4 [l,l]-4 4 [2,l]t 24 [l,2] = a 



In the above expressions e'[i, j] corresponds to ij-th el- 
ements of the matrix e'. Similarly we call the matrix 
elements of t'. 

The diagonal forms of Ml and Mr can be explained 
as follows. Due to translational invariance of the semi- 
infinite electrodes the wave amplitudes at the sites of the 
electrodes (L and R) can be written in Bloch wave form, 



A e in-PL(R) 



(21) 



Here, Pl(R) = ka, k is the wave vector and a being the 
lattice spacing of the discrete model. For the Oth site of 
the left electrode (L) we can write, 



-0ot 



J0i 



0- 



(22) 



where /3l is defined by the following energy dispersion 
relation as, 



E = e L + 2t L cos/3 L . 



(23) 



Now, the difference equation for the Oth site (the bound- 
ary site of the left electrode) is in the form, 



(E - ez,)Vo,t(4.) = *i^-i,f(4.) + *£B^i,t(40 ( 24 ) 



Thus, 
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(25) 



Using Eqs. (|22|) and (|25|) we can construct Ml which 
shows the above diagonal form. In an exactly similar 
way we get the diagonal form of Mr. 

Let us first discuss the case of up spin incidence. Con- 
sidering the whole system (left electrode-conductor-right 
electrode), the wave amplitudes at sites 0,-1 and N+l, 
N + 2 can be written in terms of reflection and transmis- 
sion amplitudes as, 



and. 



ip-H 
^04. 



V'A r +2f 
tpN+21 
IpN+lt 
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tt 
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T n e i(N+2)0R 
T n e i(N+2)p R 
T n e i(N+l)0 R 
T n e i{N+l)0 R 



(26) 



(27) 



where, represents the reflection amplitude of an up 
spin as an up spin, and p^ denotes the the reflection 
amplitude of an up spin as a down spin, r" corre- 
sponds to the transmission amplitude of an up spin with- 
out any flipping, whereas denotes the spin flip trans- 
mission amplitude. We solve the transfer matrix equa- 
tion (Eq. |H])) substituting the explicit expressions of the 
wave amplitudes from Eqs. ([2"rJ)) and (|2"T|) . 
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The transmission probabilities T-j-j- and T-^ arc defined 
by the ratio of the transmitted flux to the incident flux 
as. 



ri 



T 



t R sin (3 R 
t L sin /3 L 
t R sin /3 R 
t L sin Pi 



rtt|2 



.til 2 



(28) 



Therefore, the total transmission probability of an up 
spin becomes, 



- M 



(29) 



In a similar way, we can calculate the total transmission 
probability for the case of a down spin incidence as, 



II 



T 



It 



(30) 



Based on the Landauer conductance formula @ , the con- 
ductance g aa > through the interferometer can be calcu- 
lated. At much low temperatures and bias voltage, it can 
be expressed in the from, 



e 2 

T 

h 



(31) 



The spin dependent current flowing through the inter- 
ferometric geometry can be determined from the expres- 
sion @, 



(fs-fD)T aa ,(E)dE (32) 



studied in various aspects of interferometer-to-electrode 
coupling strength, AB flux <j>, external magnetic field 
and spin flipping. Here, we assume that the two non- 
magnetic (NM) electrodes are identical in nature. For 
our illustrative purposes, let us first mention the values 
of the different parameters those are considered for the 
numerical calculations. The on-site energies (eo) in the 
interferometer are chosen to be 0. Magnitudes of all the 
local magnetic moments (h n ), associated with the atomic 
sites of the interferometer, are fixed at 0.5 and hence- 
forth we call it simply as h. The hopping integral be- 
tween the nearest-neighbor sites of the interferometer is 
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where, fs(D) = f (E — Hs(D)) gives the Fermi distribu- 
tion function of the two electrodes with the electrochem- 
ical potential fis(D) = Ep ± eV/2. 

Due to spin flip scattering individual spin currents (1^ 
and I^i) are no longer constant quantities, whereas the 
total spin is conserved. However, we have defined spin 
currents by introducing the effect of spin-flipping. They 
are expressed as: J t = I n + I n and I± = J u + 
Non conservation of spin inside the magnetic quantum 
interferometer gives rise to a torque known as spin flip 
induced spin torque [35j . We can define this spin torque 
through the relation, 



TfUp 



d{S) 
dt 



n 



H,S 



(33) 



where, S = \a. The spin torque may produce an angular 
displacement in the orientations of the local magnetic 
moments associated with each magnetic site. 



III. NUMERICAL RESULTS AND DISCUSSION 

Spin dependent transport properties through the mag- 
netic conductor having interferometric geometry are 



FIG. 3: (Color online). Typical conductance g as a function of 
the magnetic flux <f>, where the green and magenta curves cor- 
respond to the weak- and strong-coupling limits, respectively. 
Other parameters are as follows: E = and #3 = 65 = it/3. 

set at t = 3, while, for the NM electrodes it is chosen as 
ti = t R = 4. The site energies (£l(a)) of all the sites 
in the electrodes are put to 0. The azimuthal angles tp n 
for all n are fixed at zero and also the equilibrium Fermi 
energy Ep of the conductor is set at 0. Here, we choose 
the units c = e = h = 1 for the sake of simplicity. 

Throughout the analysis, all the essential features of 
spin transport are studied for the two distinct regimes, 
depending on the strength of coupling of the interferom- 
eter to the NM electrodes. 

Case 1: Weak-coupling limit. 

This regime is typically defined by the condition 
thD(DR) << t- Here, we choose the values of the hopping 
parameters as, tip, — tp> R = 0.5. 

Case 2: Strong-coupling limit. 

This limit is described by the condition tpp>(p) R ) ~ t. In 
this regime, we set the values of the hopping strengths 
as, tpp = t DR = 2.5. 
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A. Variation of conductance with magnetic flux <f> 

As representative examples, first in Fig. [3] we plot the 
variation of and for the interferometer as a func- 
tion of magnetic flux <j>. The results are computed for 
the injecting electron energy E = 0, where the green 
and magenta curves correspond to the weak- and strong- 
coupling limits, respectively. The direction of local mag- 
netization of the magnetic atoms (labeled as 3 and 5 in 
Fig. [1} in the interferometric arms are chosen to be ori- 
ented at an angle 7r/3 with respect to the preferred +Z 
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FIG. 4: (Color online). Typical conductance g as a function 
of 6 (83 = 85 =9), where the green and magenta curves cor- 
respond to the weak- and strong-coupling limits, respectively. 
Other parameters are as follows: E = and <f> — <j)o /4. 

direction. All the other moments are aligned along +Z 
direction. Figure [3] shows that both and g\i vary 
periodically with <f> showing O flux-quantum periodicity. 
For a symmetrically connected interferometer i.e., hav- 
ing identical configuration in the upper and lower arms 
(#3 = #5), the transmission probability drops exactly to 
zero at <f> = 4>o/2 (= 0.5 in our chosen unit) and it can 
be shown very easily by simple mathematical calculation 
as follows. 

For a symmetrically connected interferometer, the 
wave functions passing through the upper and lower arms 
of the interferometer are given by, 

if J A.Sr 

ipt = ip e 71 
ip 2 = tpoe 72 



(34) 



where, ji and 72 are used to indicate the two different 
paths of electron propagation along the two arms of the 
interferometer, ipo denotes the wave function in absence 



of magnetic flux <fi and it is same for both upper and lower 
arms as the interferometer is symmetrically coupled to 
the electrodes. A is the vector potential associated with 
the magnetic field B by the relation B = V x A. Hence 
the probability amplitude of finding the electron passing 
through the interferometer can be calculated as, 



^2 1 2 



2|Vo| 2 + 2|Vo| 2 cos 



2lTC 



(35) 



where, <f> — § A.dr = J J B.ds is the flux enclosed by the 
interferometer . 

Here, it is clearly observed from Eq. (|3"5j) that at 
<fi = 4>o/2, the transmission probability of an electron ex- 
actly drops to zero. This aspect can be utilized to design 
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FIG. 5: (Color online). Up spin conductance as a function of 
energy E for the interferometer with <f> — (f>o/2 in the limit of 
weak-coupling, (a) 63 — 8 5 = 0, (b) 83 = and #5 = ir, (c) 
83 = 7T and 85 = and (d) 83 = 85 = n. 

an XOR gate which we will describe in the forthcoming 
sub-sections. The g-E spectra shows that the typical con- 
ductance gets enhanced significantly with the increase of 
coupling strength. Beside this, it is also observed that 
the transmission probability due to spin flipping is con- 
siderably smaller than the pure spin transmission. 



B. Variation of conductance with polar angle 8 of 
the magnetic moments 

In Fig. |4] we present the variation of g^ and with re- 
spect to 9, where 8 corresponds to the angle made by the 
local magnetic moments in the interferometric arms (sites 
labeled as 3 and 5) with the preferred +Z direction. The 
orientations of local magnetic moments can be changed 
by applying an external magnetic field. All the other 
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moments are aligned along +Z direction. The results 
are calculated for the typical magnetic flux = 0o/4 
and the injecting electron energy E = 0. In this case 
g-ft shows 2-7T periodicity as a function of 9 both for the 
weak- and strong-coupling limits. But, for up and down 
orientations of the local magnetic moments in the inter- 
fcromctric arms i.e., for 9 = or any integer multiple of 
7r, no spin flip takes place, and accordingly, g^ drops to 
zero for these configurations. 

Explanation of zero transmission probability for spin 
flipping is given as follow. Spin flip occurs due to the 
presence of the term h.a in the Hamiltonian (see Eq. ©), 
a being the Pauli spin matrix with components a x , a y 
and a z for the injecting electron. The spin flipping is 
caused because of the operators a+{= a x +ia y ) and c_(= 
a x — i(Ty), respectively. For the local magnetic moments 
oriented along ± Z axes, h.a (= h x a x + h y a y + h z a z ) be- 
comes equal to h z a z . Accordingly, the Hamiltonian does 
not contain a x and a y and so as a+ and cr_ , which pro- 
vides zero flipping for up or down orientation of magnetic 
moments. For this typical configuration of the localized 
magnetic moments mentioned above (±Z axis) the spin 
flip torque vanishes which is clearly seen from Eq. (|33l) . 
On the other hand, for any other orientation (apart from 
up or down configuration) of local magnetic moments 
a non-vanishing spin transfer torque appears which can 
provide angular displacements of these moments. Due to 
these displacements, an additional contribution can oc- 
cur to the spin transmission which is neglected in our 
present study. 



conductance peaks are associated with the energy eigen- 
values of the interferometer. With the increasing number 
of magnetic atoms, comprising the interferometer, num- 
ber of energy levels associated with the interferometer 
increases, therefore more resonant peaks appear in the 
conductance spectrum. 

These features can be explained as follows. When both 
the two inputs are either low (#3 = #5 = 0) or high 
(83 = #5 = 7r), the transmission probability exactly van- 
ishes at the half flux-quantum value of </> which provides 
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C. XOR gate response 

1. Conductance- energy characteristics 

Let us now describe how such a simple geometric model 
can be implemented as an XOR gate. For the forthcom- 
ing discussion, we set the AB flux <j> at </>o/2. 

In Figs. [S] and [51 we show conductance-energy (g-E) 
characteristics for the interferometric geometry both in 
the weak- and strong-coupling limits. With the help of 
external magnetic field, the orientation of the magnetic 
moments at sites 3 and 5 (measured by the parameters #3 
and 65, respectively) can be changed. Here we will show 
that, depending on the values of #3 and #5 the interfer- 
ometric geometry exhibits XOR gate response, keeping 
all the other moments oriented along +Z axis. These 
two (03 and 85) are treated as the two inputs of the 
XOR gate. Let us first discuss the case of weak-coupling 
(Fig. [5]). When both the two inputs to the gate are zero, 
i.e., 6*3 = #5 = 0, conductance exactly drops zero (see 
Fig. EJa)). On the other hand, if any one of the two in- 
puts is high i.e., 9 3 or #5 has a non-zero value, the conduc- 
tance shows fine resonant peaks (see Figs. [5jb) and (c)) 
for some particular energy values. Finally, when both 
the inputs to the gate are high, the conductance again 
vanishes (Fig. Efd)) for the entire energy range. The 



FIG. 6: (Color online). Up spin conductance as a function of 
energy E for the interferometer with <f> — 4>o/2 in the limit of 
strong-coupling, (a) 83 = 85 — 0, (b) #3=0 and 85 = ir, (c) 
83 = ty and 85 = and (d) 83 = 85 = ir. 

zero conductance for the entire energy range. The van- 
ishing behavior of conductance at (f> = <po/2 for these 
symmetric configurations of the interferometer is clearly 
understood from our earlier discussion (sub-section III A) 
If the symmetry in the orientation of the local magnetic 
moments in the two arms of the interferometer is broken 
by applying an external magnetic field i.e., #3 ^ 9$ then 
the transmission probability becomes non-zero even at 
</> = 4>a/2. Since in this case 9 3 and #5 are chosen cither 
as or 7r, no spin flipping takes place, and therefore, 
the contribution to g-f comes only from the factor g-^. 
The contribution from the spin flipping to the conduc- 
tance spectrum will be observed for any other values of 9 
(8 = 83 = 6*5) apart from and ir. Even for these orien- 
tations of the magnetic moments, the total conductance 
vanishes for the symmetric configuration, while it shows 
a finite non-zero value for the asymmetric one. Thus, we 
can conclude that the XOR gate like response will remain 
unchanged at the typical AB flux <f> = cf>o jl for any value 
of 9. In this interferometric geometry the asymmetry can 
be quantified by the term A9 (A9 = \9 3 -8 5 \). The logi- 
cal XOR gate response is obtained for any non-zero value 
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of A9, while it is best observed for the maximum value 
of A9 (A0 max = 7r), which is presented in our case. 

Another important thing to be mentioned here is that 
for 9 = or 7r, spin flip transmission does not take 
place, and accordingly, the spin transfer torque Tfu p is 
zero which does not affect the magnetization direction, 
while for any other values of 9 apart from and n, Tfu p 
is nonzero which may cause a change in the direction 
of magnetization. But the change will be the same for 
all the magnetic sites having identical magnetic moments 
and hence XOR gate feature will not get affected. Thus 
we can conclude that the violation of spin conservation 
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FIG. 7: (Color online). Up spin current / as a function of 
applied bias voltage V for the interferometer with (f> = (f>o/2 
in the limit of weak-coupling, (a) #3 = 65 = 0, (b) #3 = and 
9 5 = tv, (c) <9 3 = ty and 9 5 = and (d) #3 = 9 5 = tt. 

due to spin flip scattering should not have any signifi- 
cant impact on XOR gate response. In short, we can say 
that the spin transmission probability becomes non-zero 
if and only if the moments embedded in the interfero- 
metric arms are oriented asymmetrically. Our numerical 
results clearly justify the XOR gate response. 

In the same footing, here we also present the 
conductance-energy characteristics for the strong- 
coupling limit. The results are shown in Fig. [6] All 
the basic features in the four different choices of the two 
inputs (#3 and #5) arc exactly similar to those as pre- 
sented in Fig. [51 apart from the broadening of conduc- 
tance peaks. The contribution to the broadening comes 
from the broadening of the energy levels of the interfer- 
ometer in this strong-coupling limit. It provides a sig- 
nificant effect in the study of current-voltage (I-V) char- 
acteristics which we will describe in the following sub- 
section. It is to be noted that the conductance-energy 
spectrum for down spin is exactly mirror symmetric to 
the spectrum observed for an up spin, and accordingly, 



we do not plot the results further for down spin. 



2. Current-voltage characteristics 

All the basic features of spin dependent transport ob- 
tained from conductance versus energy spectra can be ex- 
plained in a better way through the current-voltage {I-V) 
characteristics. The current across the quantum interfer- 
ometer is computed by integrating over the transmission 
curve according to Eq. (|32p , where the transmission prob- 
ability varies exactly similar to that of the conductance 
spectrum, since we get the relation g = T from the Lan- 
dauer conductance formula (Eq. (|3II) ) with e = h = 1 in 
our present formulation. 

In Figs. [7] and |H we plot the current (I) as a func- 
tion of applied bias voltage (V) both for the symmetric 
and asymmetric orientations of the local magnetic mo- 
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FIG. 8: (Color online). Up spin current I as a function of 
applied bias voltage V for the interferometer with (f> = 0o/2 
in the limit of strong-coupling, (a) 63 = #5 = 0, (b) 63 = 
and #5 = 7r, (c) 03 — 7r and 9 5 = and (d) #3 = 9s = tt. 

ments in the upper and lower arms of the interferometer. 
The flux 4> is set at </>o/2 i.e., 0.5 in our chosen unit. 
Let us first start with the case of weak-coupling (Fig. [7]) . 
For symmetric configuration of the interferometric arms 
($3 = $5 = and 6*3 = 85 = vr), the current vanishes 
for the entire range of the bias voltage V (Figs. [3a) and 
(d)). This is due to the fact that, for these cases the 
transmission probability becomes zero for the entire en- 
ergy range as we have studied earlier. On the other hand, 
for asymmetric configuration (6*3 = 0, 85 = tt and $3 = tt, 
05 = 0) of the interferometric arms, the current is non- 
zero (Figs. [7|(b) and (c)), because of non- vanishing spin 
transmission probability. This behavior becomes much 
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more clearer from Table HI where we make a quantitative 
estimate of typical current amplitude, computed at the 
bias voltage V = 10.52. It is observed that, the current 
/ reaches the value to 0.294 only when any one of the 
two inputs are high and the other is low i.e., (#3 = and 
65 = tt or 6*3 = tt and #5 =0), while for the other two 
cases (#3 = and 8 5 = or 63 — tt and 8 5 = tt), it (I) 
gets zero. From these I-V characteristics the XOR gate 
like response is clearly visualized. In this weak-coupling 
limit, the current shows step- like behavior as a function of 
the applied bias voltage V. This is due to the presence of 

TABLE I: XOR gate behavior in the limit of weak-coupling. 
The typical current amplitude is determined at the bias volt- 
age V = 10.52. 



Input-I (63) 


Input-II (0 5 ) 


Current (I) 











TT 





0.294 





TT 


0.294 


TT 


TT 






sharp resonant peaks in the conductance spectra, as the 
current is obtained from integration method over trans- 
mission function T. With the increase in applied bias 
voltage V, the difference in chemical potentials of the 
two electrodes (/ii — ^2) increases, allowing more num- 
ber of energy levels to fall in that range, and accordingly, 
more energy channels are accessible to the injected elec- 
trons to pass through the quantum interferometer from 

TABLE II: XOR gate behavior in the limit of strong-coupling. 
The typical current amplitude is determined at the bias volt- 
age V = 10.52. 



Input-I {63) 


Input-II (0 5 ) 


Current (I) 











TT 





0.785 





TT 


0.785 


TT 


TT 






the source to drain. Incorporation of a single discrete en- 
ergy level i.e., a discrete quantized conduction channel, 
between the range (fix — ^2) provides a jump in the I-V 
characteristics. 

The case of strong-coupling is depicted in Fig. [5] A 
quantitative estimate of the current is given in Table HH 
where the typical current is measured for the bias volt- 
age V = 10.52. In this strong-coupling limit, all the basic 
features are exactly similar to that given in Table U only 
the magnitude of the output current gets enhanced to a 
value of 0.785. The non-zero currents show continuum 
behavior with respect to the change in bias voltage. As 
the sharp and discrete feature of the conductance peaks is 



lost in this strong-coupling limit, acquiring some broad- 
ening, the current changes continuously providing a much 
larger amplitude. Therefore, tuning the strength of the 
interferometer-to-electrode coupling, current can be en- 
hanced significantly keeping the bias voltage constant. 
Exactly similar kind of behavior can also be observed for 
down spin current, and accordingly, we do not show the 
results here. 



IV. CONCLUDING REMARKS 

To summarize, in the present work we have explored 
spin dependent transport through interferometric geome- 
try, penetrated by a magnetic flux <fi using transfer matrix 
formalism. A simple tight-binding framework has been 
adopted to illustrate the system, where the interferome- 
ter comprised of magnetic atomic sites is sandwiched be- 
tween two non-magnetic electrodes, namely, source and 
drain. We have calculated numerically the spin depen- 
dent transmission probability including the effect of spin- 
flip. Our numerical calculations describe conductance- 
energy and current-voltage characteristics as functions of 
the interferometer-to-electrode coupling strength, mag- 
netic flux and the orientation of the local magnetic mo- 
ments associated with each atom placed in interferomet- 
ric arms. 

First, we have observed the variation of conductance 
incorporating the effect of spin flip, as a function of mag- 
netic flux cj> showing </>o periodicity. It is noticed that, 
at half flux-quantum value of (j>, transmission probability 
drops to zero for a symmetrically connected interferome- 
ter. Next, we have studied the variation of conductance 
with 9 (considering 63 = 6§ = 6) for symmetric configu- 
ration which shows that spin flipping is blocked for 6 = 
or any integer multiple of tt. Finally, we have obtained 
XOR gate like response in g-E and I-V characteristics 
depending on the orientations of the local magnetic mo- 
ments in the upper and lower arms of the interferometer 
at <f> = 4>o/2. 

Throughout our work, we have addressed all the es- 
sential features of XOR gate operation considering an 
interferometer with total 6 atomic sites. Among them 
4 are placed at the corners to form a ring like structure 
and the rest two are coupled to the electrodes directly. In 
our model calculations, this typical number (6) is chosen 
only for the sake of simplicity. Though the results pre- 
sented here change numerically with the ring size, but 
all the basic features remain exactly invariant. The main 
point of concern is that, whether the moments in the 
upper and lower interferometric arms are oriented sym- 
metrically or not. The local magnetization direction can 
be changed by a rotation of the exchange field on the 
magnetic sites [n], HH . Change of the local moment ori- 
entations at sites 1, 2, 4 and 6 in our present geometric 
model does not make any difference to the physical fea- 
tures of the results shown above. To be more specific, 
it is important to note that, in real situation the exper- 
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imcntally achievable rings have typical diameters within 
the range 0.4-0.6 (im. In such a small ring, unrealisti- 
cally very high magnetic fields are required to produce 
a quantum flux. To overcome this situation, Hod et al. 
have studied extensively and proposed how to construct 
nanometer scale devices, based on Aharonov-Bohm inter- 
fcromctry, those can be operated in moderate magnetic 
fields j37|. 

In this work, we have calculated all the results by ig- 
noring the effects of temperature, electron-electron cor- 
relation, electron-phonon interaction, disorder, etc. Here 
we fix the temperature at OK, but the basic features 
will not change significantly even in non-zero finite (low) 
temperature region as long as thermal energy (fcsT) is 
less than the average energy spacing of the energy lev- 
els of the quantum interferometer. Over the last few 
years a lot of efforts are made to incorporate the effect 
of electron-electron correlation in the study of spin de- 
pendent transport. Electronic correlation may cause de- 
coherence among the waves passing through the interfer- 
ometric arms. But, at low temperatures the decoherence 
produced by electron-electron correlation can be limited 
and in a very recent work Montambaux et al. [38[ have 
justified it by studying electron transport for some arrays 
of connected mesoscopic metallic rings in presence of elec- 



tronic correlation. The presence of electron-phonon inter- 
action in Aharonov-Bohm interferometers provides phase 
shifts of the conducting electrons and due to this dephas- 
ing process electron transport through an AB interferom- 
eter becomes highly sensitive to the AB flux p with the 
increase of electron-phonon coupling strength [39( . In the 
present work, we have addressed our results considering 
the site energies of all the atomic sites of the interferome- 
ter are identical i.e., we have treated the ordered system. 
But in real case, the presence of impurities will destroy 
phase coherence significantly which affects the transport 
properties. In this model it is also assumed that the two 
side-attached non-magnetic electrodes have negligible re- 
sistance. At the end, we would like to mention that we 
need further study in such systems by incorporating all 
these 1 effects. 

At the end, here we have designed a spin XOR gate 
using a quantum interferometer, based on the effect of 
quantum interference, which is a classical logic gate. On 
the other hand, quantum logic gates using ring geome- 
tries have already been prop osed earlier which can be 
available in the reference [40j |. 

All these predicted results may be utilized in designing 
tailor made spintronic circuits. 



[1] P. Ball, Nature (London) 404, 918 (2000). 

[2] S. Datta, Electronic transport in mesoscopic systems, 

Cambridge University Press, Cambridge (1997). 
[3] S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. 

M. Daughton, S. von Molnar, M. L. Roukes, A. Y. 

Chtchelkanova, and D. M. Treger, Science 294, 1488 

(2001). 

[4] J. Stohr and H. C. Siegmann, Magnetism- From Funda- 
mental to Nanoscale Dynamics, Springer (2006). 

[5] J. Daughton, J. Appl. Phys. 81, 3758 (1997). 

[6] A. A. Shokri and A. Saffarzadeh, Eur. Phys. J. B 42, 187 
(2004). 

[7] A. A. Shokri and A. Saffarzadeh, J. Phys.: Condens. 

Matter 16, 4455 (2004). 
[8] M. Mardaani and K. Esfarjani, Physica E 25, 119 (2004). 
[9] M. Mardaani, A. A. Shokri, and K. Esfarjani, Physica E 

28, 150 (2005). 

[10] A. A. Shokri, M. Mardaani, and K. Esfarjani, Physica E 
27, 325 (2005). 

[11] A. A. Shokri and M. Mardaani, Solid State Commun. 
137, 53 (2006). 

[12] U. Gerland, J. von Delft, T. A. Costi, and Y. Oreg, Phys. 

Rev. Lett. 84, 3710 (2000). 
[13] D. Goldnaber-Gordan, H. Shtrikman, D. Mahalu, D. 

Abosch-Magder, U. Meirav, and M. A. Kastner, Nature 

(London) 391, 156 (1998). 
[14] S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwen- 

hoven, Science 281, 5 (1998). 
[15] H. Yu and J.-Q. Liang, Phys. Rev. B 72, 075351 (2005). 
[16] J. Konig and Y. Gefen, Phys. Rev. B 65, 045316 (2002). 
[17] F. Chi and S.-S. Li, J. Appl. Phys. 100, 113703 (2006). 
[18] Y. Imry, Introduction to Mesoscopic Physics, Oxford Uni- 



versity Press, Oxford (2002). 
[19] S. Bellucci and P. Onorato, Phys. Rev. B 78, 235312 
(2008). 

[20] S. Bellucci and P. Onorato, J. Phys.: Condens. Matter 

19, 395020 (2007). 
[21] P. Foldi, O. Kalman, M. G. Benedict, and F. M. Peeters, 

Nano Lett. 8, 2556 (2008). 
[22] P. Foldi, O. Kalman, and F. M. Peeters, Phys. Rev. B 

80, 125324 (2009). 
[23] B. Molnar, P. Vasilopoulos, and F. M. Peeters, Phys. 

Rev. B 72, 075330 (2005). 
[24] Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 
[25] G. Kirczenow, Phys. Rev. B 72, 075351 (2005). 
[26] P. A. Lee and D. S. Fisher, Phys. Rev. Lett. 47, 882 

(1981). 

[27] Q. Sun, J. Wang, and H. Guo, Phys. Rev. B 71, 165310 
(2005). 

[28] S. K. Maiti, Solid State Commun. 149, 1623 (2009). 

[29] J. Heinrichs, J. Phys.: Condens. Matter 12, 5565 (2000). 

[30] S. Chattopadhyah and A. Chakrabarti, J. Phys.: Con- 
dens. Matter 16, 313 (2004). 

[31] R. Landauer, IBM J. Res. Dev. 32, 306 (1988). 

[32] M. Biittiker, IBM J. Res. Dev. 32, 317 (1988). 

[33] A. Ghosh and S. N. Karmakar, Phys. Rev. B 58, 2586 
(1998). 

[34] C. E. T. Goncalves da Silva and B. Koiller, Solid State 

Commun. 40, 215 (1981). 
[35] Z.-G. Zhu, G. Su, B. Jin, and Q.-R. Zheng, Phys. Lett. 

A 306, 249 (2003). 
[36] X. Qian, X. Gu, and Y. Ji, Rev. Sci. Instrum. 81, 106106 

(2010). 

[37] O. Hod, R. Baer, and E. Rabani, J. Am. Chem. Soc. 127, 



11 



1648 (2005). 266803 (2006). 

[38] C. Texier, P . Delplace, and G. Montambaux, [40] P. Foldi, B. Molnar, M. G. Benedict, and F. M. Peeters, 

larXiv:0907.3133V l. Phys. Rev. B 71, 033309 (2005). 

[39] O. Hod, R. Baer, and E. Rabani, Phys. Rev. Lett. 97, 



